Preprint typeset in JHEP style - HYPER VERSION 



DFTT 7/2002 



bep-ph/0202237| 



The contribution of the four-parton final states to 
7*7* — ► hadrons 



(N 
O 

o 

(N ; Vittorio Del Duca 

Istituto Nazionale di Fisica Nucleare, Sez. di Torino 
via P. Giuria, 1 - 10125 Torino, Italy 



E-mail: delduca@to.infn.it 



Fabio Maltoni 

Dept. of Physics, University of Illinois at Urbana- Champaign 
Urbana, IL 61801, USA 



(N 
> 

[~~. , E-mail: pialtoni@uiuc . edu 

m 

CN ' Zoltan Trocsanyi* 

(N : 

C^} University of Debrecen and 

Institute of Nuclear Research of the Hungarian Academy of Sciences 
H-4001 Debrecen, PO Box 51, Hungary 
' E-mail: zoltan@zorro.atomki.hu 

a: * a 

Abstract: In the analysis of the total cross section for the 7*7* — ► hadrons process, 
we include the four parton final states, which are part of the 0(ot 2 s ) corrections. The 
four-parton final states contain the diagrams with gluon exchange in the crossed channel, 
^ ■ which constitute the leading order of the BFKL resummation. We show that the diagrams 

with gluon exchange in the crossed channel play an important role in the large Y region, 
however their contribution to the cross section must be evaluated exactly. In fact, the 
high-energy limit, which constitutes the kinematic framework of the BFKL resummation, 
is not sufficiently accurate at LEP2 energies. The inclusion of the diagrams with gluon 
exchange in the crossed channel reduces the discrepancy between the theory and the LEP2 
data collected by the L3 Collaboration, but the data still lie above the theory, even allowing 
for a large scale uncertainty in the theory. Thus, in order to describe accurately the data 
for 7*7* — > hadrons in the large Y region, corrections of an order higher than 0{a 2 s ) seem 
to be necessary. 
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1. Introduction 

Strong interaction processes, characterised by a large kinematic scale, are described in per- 
turbative QCD by a fixed-order expansion of the parton cross section in a s , complemented, 
if the scattering process is initiated by strong interacting partons, with the Altarelli-Parisi 
evolution of the parton densities. However, in kinematic regions characterised by two very 
different hard scales, a fixed-order expansion might not suffice: large logarithms of the 
ratio of the kinematic scales appear, which may have to be resummed. In processes where 
the centre-of-mass energy s is much larger than the typical momentum transfer t, the sub- 
process which features gluon exchange in the crossed channel, and that usually appears 
at 0(a|), tends to dominate over the other sub-processes. That sub-process constitutes 
the leading-order term of the BFKL equation, which is an equation for the Green's func- 
tion of gluon exchanged in the crossed channel. The BFKL equation [l], ^, |3| resums the 
logarithms of type bn(s/|t|). 

Over the last decade, several observables, like the scaling violations of the i*2 structure 
function j|, ||, forward-jet production in DIS f|— [11], dijet production at large rapidity 



intervals [f[^] — fT9|, and 7*7* — > hadrons in e + e collisions [^] — p6| have been proposed 
in the literature as candidates for the detection of the BFKL evolution, and have been 
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measured and analysed as functions of observables, which aim to single out large logarithms 
of type ln(s/|i|). However, from a phenomenological point of view, in order to claim 
detection of BFKL gluon radiation in a given process in an unambiguous way, we must 
rule out any explanation of that process in terms of a fixed order expansion, or in terms of a 
different resummation. Thus, in order to make a sound BFKL analysis, we must ascertain 
first of all if: 

• the sub-process with gluon exchange in the crossed channel, i.e. the leading order of 
the BFKL resummation, dominates over all the other sub-processes; 

• the acceptance cuts of the experiment under consideration allow us to reach the 
kinematic region of the high-energy limit, where the approximations needed for a 
BFKL analysis are valid. 

The goal of this paper is to analyse whether the two conditions above are fulfilled in 
the 7*7* — > hadrons process at LEP2. Namely, we consider 

7* + 7* — ► hadrons, (1-1) 

in e + e~ collisions at photon virtualities q 2 = —Q 2 < 0, and for large centre-of-mass energies 
squared W 2 = (q\ + q^) 2 , with % being the momenta of the photons. In practice we can 
realise the scattering (1.1) in the process 



e + + e — > e + + e + hadrons, (1-2) 

of which Eq. ([O]) constitutes a subset. Other contributions to the process in Eq. ( |i~^ ) 
are, for example, those in which the incoming e + e~ pair annihilates into a photon or a Z 
boson, eventually producing the hadrons and a lepton pair, or those in which one (or both) 
of the two photons is replaced by a Z boson. However, it is not difficult to devise a set of 
cuts such that the multiperipheral process 

e + + e~ — ► e + + e~ + 7* + 7* 

(1.3) 

hadrons, 



gives the only non-negligible contribution to the process in Eq. (1.2). One can tag both 
of the outgoing leptons, and retain only those events (thus termed double-tag events) in 
which the scattering angles of the leptons are small: in such a way, the contamination 
due to annihilation processes is safely negligible. Furthermore, small-angle tagging also 
guarantees that the photon virtualities are never too large (at LEP2, one typically measures 
Q 2 ~ 10 GeV 2 ); therefore, the contributions from processes in which a photon is replaced 
by a Z boson are also negligible. Thus, it is not difficult to extract the cross section of 
the process 7*7* — * hadrons from the data relevant to the process in Eq. ( |1.2j ). Double-tag 
events have in fact been studied by the CERN L3 and OPAL Collaborations, at various 



e+e" centre-of-mass energies (y/s = 91 and 183 GeV [§0|, and 189-202 GeV f|T], ||]). 



The process ( |1.1| ) has been analysed at leading order 2J] and at next-to- leading or 



der |2^] (NLO) in a s . However, at the high end of the W spectrum, the NLO prediction 
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does not suffice to describe the data. In this paper, we consider the part of the next-to- 
next-to-leading order (NNLO) corrections which yields the dominant contribution to the 
total cross section in the large- W region. As we shall argue below, that contribution comes 
from four quark final states. Since four parton final states do not yield per se a finite contri- 
bution to the total cross section, we consider a subset of them, those with gluon exchange 
in the crossed channel, which are finite, and argue that they yield the most important 
contribution to the total cross section in the large-l^ region. 

The paper is organised as follows: in Section || we set the theoretical framework; 
in Section ||] we consider the four-quark contribution to the 0{a 2 s ) corrections, exactly 
and in the high-energy limit; analysing various theoretical predictions for rapidity and W 
distributions in Section ^, we substantiate our claim that the most important contribution 
to the total cross section in the large-VF region comes from four-quark production with 
gluon exchange in the crossed channel, and also show that the high-energy limit is not 
sufficiently accurate at LEP2 energies; in Section |5| we present our phenomenological results, 
by comparing our predictions to the L3 data |2^] (we shall not perform a comparison with 
the OPAL data, which have much poorer statistics); in Section || we draw our conclusions. 



2. The theoretical framework 

At leading order, the multiperipheral process ( |0| ) is modelled by the partonic subprocess 
— > qq : depicted in Fig. |l|(a), which has been computed for massless and massive 
final-state fermions psfl . In Ref. the 0(a s ) QCD corrections to the process ( |L~3| ) 
were computed for final-state massless quarks (sample diagrams are given in Fig. |j](b)-(c)). 
However, the NLO analysis of Ref. [^] yields a cross section behaving as 

<7 7 . r ~ 1/W 2 , (2.1) 

modulo logarithmic corrections. Thus, it is only propaedeutic to the BFKL resummation, 
whose leading-order term is based on the exchange of a gluon in the crossed channel, 
which appears only at 0(a 2 ) (Fig. |l](d)). The BFKL resummation then builds up gluon 
emission along the gluon exchanged in the crossed channel (the first rungs of the ladder 
are represented in Fig. [j](e)-(f)). The diagrams represented by Fig. |l](d) are expected to 
yield a cross section which, away from the threshold and the kinematic limit, is weakly 
dependent on W. The additional gluon emissions build up the logarithmic corrections 
which the BFKL theory resums, so that the full cross section is expected to behave as 

oo oo oo 

(T 7*7* ~ '^2 a 0j»s + aia 2 s '^2(a s Ly + a 2 a| ^ a s (a s Ly H , (2.2) 

j=0 j=0 j=0 

where L = \og(W 2 / fjL%f) is a large logarithm, and the quantity is a mass scale squared, 
typically of the order of the crossed-channel momentum transfer and/or of the photon 



virtualities. In Eq. (2.2), the second and third sums collect the contributions which feature 
only gluon exchange in the crossed channel, the second (third) sum resumming the BFKL 
(next-to-)leading logarithmic corrections; the a\, a 2 coefficients behave like l//^ w . The 
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b) 



c) 






Figure 1: Sample of diagrams contributing to the production of hadrons in the collision of two 
off-shell photons. 



ellipses refer to logarithmic corrections beyond the next-to-leading accuracy. The first sum 
in Eq. ( |2.2D is a fixed-order expansion in a s starting at O(a s ), and collects the contributions 
which do not feature gluon exchange in the crossed channel; the aoj coefficients behave like 
1/W 2 *. Thus, it is clear that the second and third sums of Eq. ( |2,2| ) will eventually 
dominate over the first sum in the asymptotic energy region W — ► oo. The second sum of 
Eq. ( |2.2[ ) has been analysed in the region W 2 3> /U^, by computing in the high-energy limit 
the a\ coefficient in the massless and in the massive [27] case. As mentioned above, 



the aoo term in the first sum has been computed in Ref. 28] for massless and massive final 



state quarks, while the aoi term has been computed in Ref. [29|, for massless final-state 
quarks. In the next paragraph, we shall illustrate that at present a calculation of the ao2 
term is unfeasible. In this work, we compute exactly the a\ coefficient in the massless limit, 
and add it to the aoo an d aoi terms. However, we do not perform the resummation, i.e. 
we consider only the j = term in the second sum of Eq. ( |2.2[ ). 

In examining the radiative corrections to the process Ql-3|), we first note that at present 



a calculation of the full 0(a s ) corrections to the total and to the inclusive jet and dijet 
cross sections is unfeasible, since it would require a computation of two-loop amplitudes, 



*The ao2 coefficient may feature terms which behave like l/(WfJ,v/) and arise from the interference 
between diagrams with gluon exchange in the crossed channel and diagrams with quark exchange in the 
crossed channel. These interference terms have been analysed in Section [i] (see the discussion of Fig. ^ on 
page 111]). 
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b) 




d) 



Figure 2: Amplitude for e + e" 



e) 



e + e qqQQ scattering. Figures (a), (6) and (c) 



correspond to the functions g a , gb and g c in Eq. (3.3); figures (d), (e) and (/) contribute to the 
function /. 



including double box diagrams with two off-shell legs which are not known. A calculation 
of the 0(a 2 s ) corrections to the inclusive three-jet rate is feasible but difficult, and given 
the very limited experimental statistics available it would have at present only an academic 
value. Therefore, we shall limit ourselves to the analysis of the contribution of the four 
parton final states to the 0{a 2 s ) corrections. In particular, the final states can be made of 
either two quark pairs or a quark pair and two gluons. In Fig. |2] the diagrams with two 
final-state quark pairs are represented. 

The treatment of the four parton final states poses some additional problems, be- 
cause as far as the total, the inclusive jet and the dijet cross sections are concerned, only 
the contribution of the four-quark diagrams with gluon exchange in the crossed channel 
(Fig. ^(a)-(c)) is infrared finite. Henceforth, we shall term the diagrams with gluon ex- 
change in the crossed channel the g class. The diagrams of the g class are by themselves 
gauge invariant and as described in the former paragraph, they are expected to yield a 
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cross section which is only logarithmically dependent on W 2 . Thus, as W 2 grows, they 
are expected to dominate over the diagrams with quark exchange in the crossed channel. 
The diagrams with two final-state quark pairs and quark exchange in the crossed channel 
(Fig. 0(d)- (f)) and the diagrams with a quark pair and two gluons in the final state, present 
infrared divergences when one or both the final-state gluons become collinear to the parent 
quarks or between themselves, when two or three final-state quarks become collinear, or 
when the gluon emitting the quark pair becomes soft. We shall term the diagrams with a 
quark pair and two gluons in the final state, and the diagrams with two final-state quark 
pairs and quark exchange in the crossed channel the / class. The diagrams of the / class 
would all be part of the proper NNLO corrections to the process 7*7* — ► qq, and after 
cancellation of their infrared divergences by the corresponding virtual terms their contri- 
bution to the cross section is expected to have a dependence on W 2 like in Eq. ( |2.1| ), up 
to logarithmic corrections. They would contribute the CI02 coefficient in Eq. (2.2). 



3. The four-quark contribution to the 0(ct 2 ) corrections 

When calculating the four-quark contribution to 7*7* — ► hadrons, 

e + e~ — > e + e~7*7* — ► e + e~qqQQ, (3-1) 

we assume that all the produced quarks are massless and that the virtualities Q 2 2 are low 
enough that contributions of virtual W, Z bosons can be neglected. Indeed, for the four- 
quark contributions at 200 GeV centre-of-mass energy (W) ~ 40GeV and (Qf) ~ 13GeV 2 . 
We have re-computed the tree amplitudes for four-quark production |K| using the spinor 
products of Appendix The contributing diagrams are shown in Fig. ||, featuring gluon 
(Fig. ||(a) — (c)) or quark (Fig. ||(d) — (/)) exchange in the crossed channel. For two quark 
pairs of different flavour, we have 

^ 8 (1 9 ,2 9 -;3,,4,-,5,-,,6^;7q,8q) =4e 4 g 2 TfeTfa A s (l, 2; 3, 4, 5, 6; 7, 8) , (3.2) 

with {1, 2} and {7, 8} the quark pairs, and {3, 4} and {5, 6} the lepton pairs, and where the 
colour-stripped sub-amplitude A$ depends on the momenta and helicities of the external 
particles. By convention, all particles are taken as outgoing, thus an incoming fermion of 
a given helicity is represented by an outgoing antifermion of the opposite helicity. A% can 
be divided into the functions as, bs and eg, 

A 8 (l q , 2 q ; 3 e , 5|>, 6 e ; 7q, 8q) = Q fq Qf Q a 8 (l, 2; 3, 4, 5, 6; 7, 8) 

+Q%h(l, 2; 3, 4, 5, 6; 7, 8) + Q%c 8 (l, 2; 3, 4, 5, 6; 7, 8) , (3.3) 

with Qf q(Q) the electric charge fraction of the quark q(Q) of flavour / g (Q). The calculation 
of the functions a§, bs and c 8 as well as the one of the related production rate is detailed 
in Appendix |B[ 
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3.1 The high-energy limit 

In Section [2], we claimed that when the squared hadronic energy W 2 is much larger than the 
typical momentum transfer /j 2 ^, and radiative corrections to the gluon exchanged between 
two quark pairs are considered, the large logarithms of type ln(W 2 / fj^) which ensue can 
be resummed through the BFKL ladder. In fact, the resummation of the BFKL ladder 
requires more restrictive kinematics, where the rapidities of quarks which do not belong to 
the same quark pair are strongly ordered, 

771 ~ r/ 2 > 777 ~ r/ 8 or 771 ~ 772 <C 777 ~ r] 8 . (3.4) 

Eq. ( |3.4D defines the high-energy limit for 7*7* — ► hadrons. When the strong rapidity 
ordering ( |3.4j) occurs, the diagrams with gluon exchange in the crossed channel yield the 
dominant contribution^ , and in the amplitude the functions 6s and c 8 can be neglected. 
Thus the two final-state quark pairs can be treated as non-interacting, and in the squared 
amplitude we can take the quark flavours as always distinct, 

i J2 \Mh 2; 3, 4, 5, 6; 7, 8)| 2 = \{Q\n u + Q 2 d n d ) 2 \a 8 (l, 2; 3, 4, 5, 6; 7, 8)| 2 , (3.5) 

fqjQ 



with Q u = 2/3, Qd = —1/3 and n u u\ the number of up(down)-type quarks. In Eq. (3.5) the 
factor of 1/2 appears in order to avoid double counting, since the amplitudes are symmetric 
with respect to the interchange of the two quark lines. Thus for a fixed lepton-helicity 
configuration, e.g (3^T, 4J, 5^~, the production rate is 

<M37,4+ 5f„6J,) (3.6) 
= ^dP 6 16 (N 2 - 1) (4vra em ) 4 (47ra s ) 2 i(Q 2 n u + Q 2 d n d ) 2 

a 8 (l-,2 + ;3-,4 + ,5 + ,6-;7-,8 + )| 2 + |a 8 (2-,l + ;3-,4 + ,5 + ,6-;7-,8 + )| 2 
+ |a 8 (l-,2+;3-,4 + ,5 + ,6-;8-,7 + )| 2 + |a 8 (2-,l + ;3-,4 + ,5 + ,6-;8-,7+)| 



where the six-particle phase space is given in Eq. ( B.15 ). The other lepton-helicity config- 
urations are simply obtained by exchanging the labels 3 and 4 and/or 5 and 6 in Eq. (fTfj) 
(see Appendix |B|). The unpolarised rate is given by averaging over the rates for the four 
lepton-helicity configurations. 

Since in the high-energy limit the two final-state quark pairs behave effectively as 
if they were two independent scattering centres, the amplitude fl3.3|) with the functions 
6 8 and c 8 set to zero is expected to factorise into two high-energy coefficient functions, 
usually termed impact factors, for the process eg* — > eqq, where g* is the off-shell gluon 
which is exchanged in the crossed channel. In the high-energy limit, the amplitude ( p.5| ) 
can then be used to derive such impact factors. However, it is easier to invoke high-energy 



^Note that the reverse is not true. In fact the diagrams with gluon exchange in the crossed channel may 
dominate over the diagrams with quark exchange well before the high-energy limit is realised. This issue is 
discussed in Section \% 
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factorisation and to derive the impact factor for eg* — > eqq from a simpler process, e.g. 
from the scattering amplitudes eg — > qqg or eQ — ► qqQ which appear typically in DIS 
processes. Then one can use two such impact factors, one in the forward and one in the 
backward kinematics, connected by a gluon exchanged in the crossed channel, in order to 
obtain the amplitude for e + e~ — * qqQQ in the high-energy limit. We denote the impact 
factor for the eg* — > eqq process in the forward (backward) kinematics, evaluated in the 
7*7* centre-of-mass frame, as Vf^(pfpe,p q ,pq) (remember that all momenta are outgoing, 
hence the dependence on momenta of two particles and two antiparticles) , and derive it in 
Appendix [C] (see Eq. flC.16|) ). Then the amplitude (3.2) factorises as 



•4-8 2 5 ; 3 e , 5j/ , 6 e r ;7q,8q) 
= 2s 



l r 



9se 2 Q fq Tfc V2V f (A f 3 e , l q , 2 9 )\ - [g s e 2 Q fQ If rh V2V b (6 e ,; 5 £ -„ 7 Q , 8 Q )\ (3.7) 



with t = q 2 , where q = — X)i=iPi i s the momentum transfer. The two impact factors for 
eg* — > eqq can be extracted through re-labelling from Eqs. ( |C.16| )-( C.19| ). Each of those 
impact factors can be decomposed further into a lepton current and an impact factor for 
1*9* -» qq (see Eqs. (|O20|) -(|C^) ). 

In computing the square of the amplitude, we must sum over helicity, colour and flavour 
of the quarks, however in this case the flavour sum is trivial since the two impact factors 
do not interfere and we can treat the quark flavours as always distinct. The production 
rate is 

da(37,4+,5+,6^) = ^dV 6 16 (N 2 C - l) {^a Gm )\AKa s ) 2 {Q 2 u n u + Q 2 n d f 



'7*7* 
t 2 



\Vf{4+;3J,l-,2+)\* + \V f (4p,3J,l+,2i 
x (|y b (5t;6 £ -;,7 Q ,8t)| 2 + |^(5t;6,-,7+85)| 2 ) , (3. 



which constitutes the high-energy factorisation of Eq. (|3.6| ). Each of the two rapidity 
orderings of Eq. (3^) yield the same contribution to Eq. ( |3.8|) . Thus we have included 
them by taking only the first of the two and deleting the double counting factor 1/2. As in 



Sect. 2.3 of Ref. [^9|, the phase space ( B.15 ) for the e + e qqQQ final state can be factorised 
into hadronic and leptonic phase spaces, 



dV G = dT(p 3 ,p 5 )d'P4(pi,P2,P7,P8;h +k 2 ) 



(3.9) 



with ki = pi — ps and k 2 = pe — p§ the momenta of the virtual photons (here we have 
inverted the direction of p^ and p% in order to have them incoming) , and 



dr 

dV 4 



d 3 p 3 d 3 p 5 
(2vr) 3 2p0 (2vr) 3 2p° 

II T^Askfl ( 2?F ) 4 * 4 ( fc l + fc 2 - Pl - P2 - P7 ~ PS) , 



J- J- (2vr) 3 2pC 

=1,2,7,8 v ' ^ 

the leptonic and hadronic phase spaces, respectively. 



(3.10) 
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In the high-energy limit, momentum conservation for e + e — * e + e 7*7* — > e + e qqQQ 
implies that, in the 7*7* centre-of-mass frame^, 

k + ^ p + + p + 

k 2 -P7 +P 8 (3-11) 
= Pijl + P2± + P7± + Ps± 

where we use light-cone coordinates: p^ = p° ± p z and for the two-dimensional vector p± 
complex transverse coordinates p± = p x + ip v (see Appendix |C]) . Momentum conservation 
fl3.11 ) allows us to factorise the hadronic phase space fl3.10[) further, 



x {2tt) z 5 z ( Pl± + p 2 ± + P7± + P8±) ■ (3.12) 

The terms in round brackets in the first line are the phase spaces for the two impact factors. 
They are connected by transverse momentum conservation only. The overall factor of 2 in 
dVi comes from the Jacobian of the light-cone coordinates. Eq. ( 3.12] ) can be immediately 



generalised to the emission of a BFKL gluon ladder between the impact factors. 
Fixing 

X a = T TT = 1 %a j Xf) = — — = 1 X\) , (3.13) 

Pi +P2 P7 + Ps 

the phase space ( |3 . 1 2| ) can be re-written as 

11/ dx a d 2 pi±\ ( dx b d 2 p 71 _ 

oPa 



(4tt) 2 2A^*2 \x a (l-x a ) (2tt)V V^(l-^b) ( 27r ) 2 

(2vr) 2 5 2 (( ?a± - %± ) (3.14) 



(2vr) 2 (2tt) 5 

with q a = k% — pi — P2 and q^ = pj + ps — &2- Note that in the 7*7* centre-of-mass 
frame, the momenta of the virtual photons are (in the light-cone notation of Appendix |C|) 
k\ = (kf ,k± ;0±) and k2 = (k%, k^', 0j_), with virtualities /c 2 = kfk^[ = —Q\ and fc| = 
k^k^ = — Q 2 ,. I n the high-energy limit, kf 3> and k^ 3> J, thus the centre-of-mass 
energy is s 7 * 7 * = (k\ + /C2) 2 ~ kfk^- 

4. Theoretical predictions 

In this section we present the results obtained by considering the contribution of the four 
parton final state to the cross section for 7*7* — > hadrons. The four parton final state is 
0{a^ ra a 2 s ), however from the stand point of both the electromagnetic and the strong cor- 
rections it is a leading order calculation, thus the dependence of either the electromagnetic 
or the strong coupling on the respective scales is maximal. 



*Eq. (3.11) is valid also in e + e centre-of-mass frame by adding p$± +p$± to the right hand side of the 
third line. 
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As far as a em is concerned, we have chosen to set the scales on an event-by-event 
basis to the virtualities of the exchanged photons; hence, we replace the Thomson value 
«o — 1/137 by a em (Qf), as in Ref. [29]. This choice better describes the effective strength 
at which the electromagnetic interaction takes place. In addition, we treat independently 
the two photon legs: thus, in the formula? relevant to the cross sections, a^ m has to be 

Understood as a 2 , m (Q\)a 2 , m {Q\) ■ 

As far as a s is concerned, we define a default scale /jlq so as to match the order of 



magnitude of the inverse of the interaction range |29], 



2 _Qt + Q 2 . f Pi± + P2± + P7± + P8± 



^ = ^y^ + [ — — 2 " ' — ) • d-i) 



Scale choices other than (4.1) have been considered in Ref. p9| . The renormalisation scale 



/i entering a s is set equal to //q as a default value, and equal to fJ-o/2 or 2/io when studying 
the scale dependence of the cross section. In Eq. ( [4.1[ ), the are the transverse energies 
of the outgoing partons. Since the hard process is initiated by the two virtual photons, 
the proper frame to study its properties is the 7*7* centre-of-mass one. Therefore, when 
talking about transverse energies, whether in a total or in a jet cross section, this frame 
will be always understood. 

We evolve a s to next-to-leading log accuracy, with a s (M z ) = 0.1181 |32| (in MS at two 
loops and with five flavours, this implies = 0.2275 GeV). The choice of the two-loop 
running is due to the fact that we are going to use a full NLO calculation augmented by a 
partial 0(a 2 s ) contribution (the diagrams of the g class only). When presenting numerical 



results we use five massless flavours in the cross section formulae ( |3.&j ) and (B.9). 

In exploring the footprints of the BFKL resummation in e + e~ collisions, it is customary 
to introduce the variable Y, 

F = log ^ a (4.2) 
V Q1Q2 

where the variables yi are proportional to the light-cone momentum fraction of the virtual 
photons, 

9i + if , 2 6i 
V%= 7= = 1 1= cos — , z = l,2, (4.3) 

where Ei and Q{ are the energies and scattering angles of the outgoing electron and positron 
in the e + e~ centre-of-mass frame. For large Y, we have y\y 2 s ~ W 2 , i.e. the Y variable 
parametrises the ratio of the hadronic energy over a typical momentum transfer, thus it is 
a variable which is suitable for analyses of the BFKL type. 

In Fig. ^ on the next page, we plot the Y distribution for different production rates. 
Namely, for the total NLO cross section (dot-long-dashed line), for the total cross section 
due to diagrams of the g class only (solid line), for the four-jet cross section at y cut = 0.01 or 
0.001 (dot-short-dashed line) and for the four-jet cross section with diagrams of the g class 
only at y cnt = 0.01 (dashed line) and at y cn t = 0.001 (dotted line). Throughout this plot 
and henceforth, we use the acceptance cuts of the CERN L3 Collaboration f22|| , namely, 
the lepton energies Ei are larger than 40 GeV, the lepton tagging angles 9i are between 
30 and 66 mrad and the hadronic energy W is larger than 5 GeV. As far as the photon 
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Figure 3: Y distribution in different production rates. Namely, in the total NLO cross section 
(dot-long-dashed line), in the total cross section due to diagrams of the g class only (solid line), 
in the four-jet cross section (dot-short-dashed line) and in the four-jet cross section with diagrams 
of the g class only at y cut = 0.01 (dashed line) and at y cut = 0.001 (dotted line). The relative 
normalisation of the curves has been rescaled, in such a way that the area under each of the curves 
is the same. 



virtualities are concerned, the cuts above imply that Q\ 2 > 4 GeV 2 . For the cross section 
due to diagrams of the g class only, we use Eq. (|3.6|) . For the four-jet cross section, the four 
quark final states have been computed through the formulae of Appendix the final states 
with two quarks and two gluons have been generated with the help of MADGRAPH |}3| . 
which has been used also to check numerically the four-quark amplitudes of Appendix [B]. 
In the jet cross sections, we define the jets through a k T algorithm |]3j}] . The jet size is 
set by the y cu t variable. In Fig. ^ and |I| the normalisation of the curves is not relevant, 
since the contribution of the diagrams of the / class to the total cross section cannot be 
inferred from the four-jet cross section, due to the lack of virtual corrections. Thus the 
relative normalisation of the curves has been rescaled, in such a way that the area under 
each of the curves is the same. Note that in Fig. || the shape of the four-jet cross section 
(dot-short-dashed line) is largely independent of the chosen y cu %. In fact, there is basically 
no difference between the dot-short-dashed line at y cn t = 0.01 or 0.001. At large Y the 
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Figure 4: Rapidity distribution of one of the final-state partons, in different production rates. 
Namely, in the total LO cross section (circles), in the total cross section due to diagrams of the g 
class only (diamonds), in the four-jet cross section (circles with crosses) and in the four-jet cross 
section with diagrams of the g class only (crosses). 



four-jet cross section has a similar shape as the total NLO cross section (dot-long-dashed 
line). On the contrary, at large Y the four-jet cross section with diagrams of the g class only 
gets a larger and larger contribution as y cut goes from 0.01 (dashed line) to 0.001 (dotted 
line). In addition, the diagrams of the g class are by themselves infrared finite. Thus in the 
four-jet cross section from diagrams of the g class only, we can take the limit y cu t — * and 
obtain the total cross section (solid line). This has the most open shape at large Y, which 
hints that at large Y we should expect a substantial contribution from the diagrams of the 
g class to the total cross section at 0(a 2 s ). The interference terms between the diagrams of 
the g class and those of the / class, i.e. the terms 2[Q^Qj Q Re(a|6g) + Q f q Q'j Q R,e(alcs)] in 
Eq. (B.9), give also a finite contribution to the total cross section. We have checked that 
they yield a curve that is similar in shape to the total NLO cross section. Compared to 
the diagrams of the g class (solid line), they yield an increase of at most 10% only in the 
small-l" region, their contribution being negligible at large Y. Thus we shall neglect them 
henceforth. 

In Fig. ||, we plot the rapidity distribution of one of the final-state partons, in different 
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production rates. Namely, in the total LO cross section (circles), in the total cross section 
due to diagrams of the g class only (diamonds), in the four-jet cross section (circles with 
crosses) and in the four-jet cross section with diagrams of the g class only (crosses). Note 
that, as expected, in the total LO cross section [29|, which features two final-state partons 
with quark exchange in the crossed channel, the partons produced populate mainly the 
central rapidity region (circles). That is true also for the four-jet cross section (circles with 
crosses). Instead, in the four-jet cross section from diagrams of the g class only (crosses), 
at Vcut = 0.001, the produced quarks populate mainly the forward and backward rapidity 
regions. However, the shape of the circles with crosses curve depends on y cu t very mildly, 
while that of the crosses curve depends strongly on the chosen y cu t- To understand how 
the latter comes about, we recall that the diagrams of the g class feature two quark pairs 
separated by a gluon exchanged in the crossed channel, and therefore susceptible of being 
produced at large rapidity. In fact, the probability of finding one of them in the forward 
and backward rapidity regions grows as y cn t becomes smaller. In the limit y C ut — ► we 
obtain the total cross section from diagrams of the g class only (diamonds). That shows 
that it is more likely to produce the quarks in the forward and backward rapidity regions 
than it is to produce them in the central region^. Finally, we recall that in a full 0(a 2 s ) 
calculation of the total cross section, at present unfeasible, the diagrams of the / class 
are expected to yield a small correction to the NLO total cross section. Thus in a full 
NNLO calculation of the total cross section in the large Y region, we expect the rapidity 
distribution of one of the final-state partons to be roughly a combination of the circles and 
diamonds curves. 

Figures || and ||] show that in the large Y region, the diagrams of the g class yield an 
important contribution to the total cross section. Then it is natural to ask if, within the 
acceptance cuts of the LEP2 Collaborations, the hadronic energy W is sufficiently high to 
warrant the use of the high-energy limit ( |3.4| ) . We can answer that by comparing the exact 
contribution of the diagrams of the g class ([O]) to the high-energy limit of the squared 
matrix element integrated over the exact phase space ( |3.8| ). That comparison is shown in 
Fig. [| on the next page, where the solid line is the contribution of the diagrams of the g 
class and the dotted line is the high-energy limit of the squared matrix element. In this and 
in the following plots, the high-energy limit is obtained in the 7*7* centre-of-mass frame. 
In order to match the experimental accuracy, in the theoretical prediction we consider 
the high-energy limit as accurate only if the difference between the exact calculation and 
the high-energy limit is less than 20%, which is a conservative upper limit of the total 
experimental error in the tail of the distributions (see the last bins in Table |]). In Fig. || 
we see that their difference is less than 20 % only for Y > 7. Unfortunately, the region 
where Y > 7 is negligible at the LEP2 experiments, for the kinematic limit of Y ~ 8 
is almost reached, thus the statistics are very small. We can collect the events in the 
high-energy region by separating the forward and backward rapidity regions, which can be 
achieved by requiring that the sum of the rapidities of the two most forward momenta, 77/, 
be larger than 3 and that of the other two (backward) momenta, 775, be less than —3. The 

^The position of the peaks as well as the depletion in the central rapidity region depends also on the cut 
on Y, e.g., if we use the cut Y > 5, the peaks move to about ±3 and become more pronounced. 
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Figure 5: The total cross section as a function of Y. The solid (dashed) line represents the 
contribution of the diagrams of the g class (with a rapidity cut), the dotted (dot-dashed) line 
represents the high-energy limit of the squared matrix element (with a rapidity cut). 



corresponding cross sections are also shown in Fig. [5] for the diagrams of the g class (dashed 
line) and for the high-energy limit of the squared matrix element (dot-dashed line) . These 
curves almost coincide with the cross sections without the rapidity separation for Y > 7, 
but are distinctively smaller for Y < 7. This is consistent with the difference between the 
solid and the dotted lines above. In addition, it shows that for a realistic set-up, i.e. for the 
cuts of the L3 Collaboration, the high-energy limit ( |3.4| ) is more stringent than the limit 
W 1 S> the difference between the two being numerically negligible only for Y > 7. 
Since the high-energy limit ( |3.4| ) is the kinematic framework of the BFKL resummation, 
we conclude that if a BFKL resummation is used in the Y < 7 region, we expect the 
subleading logarithmic corrections to be sizeable. 

In Fig. [6| on the following page, we plot the cross section as a function of the hadronic 
energy W using different approximations: the solid line is the contribution of the diagrams 
of the g class (^^) , the dashed line and the dot-dashed line are the high-energy limit of the 
squared matrix element integrated over the exact phase space (|3.8| ) and the high-energy 
phase space (|3.14 ), respectively. In the last one, the limit s 7 * 7 * — > 00 is taken and the 
transverse momenta of the quarks are integrated out analytically. Thus we shall term it 
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Figure 6: The total cross section as a function of the total hadronic energy W. The solid line 
is the contribution of the diagrams of the g class, the dashed line and the dot-dashed line are the 
high-energy limit of the squared matrix element integrated over the exact phase space and the 
high-energy phase space, respectively. 



the analytic high-energy limit. However, in this case we cannot use a renormalisation scale 



like Eq. (4.1), which is defined on an event-by-event basis. Thus, for the sake of comparison, 



each cross section in Fig. |6| was obtained using the scale 

A = • (4.4) 

The three curves converge only for W > 100 GeV, which is again negligible at LEP2 for the 
poor statistics of the data. In addition, we see that for W < 100 GeV the analytic high- 
energy limit ^ (dot-dashed line) significantly overestimates the exact contribution of the 
diagrams of the g class. For instance, in the [40, 100] GeV range, considered by the CERN 
L3 Collaboration [^] (see the next section), the analytic high-energy limit prediction is 
about 60 % larger than the exact four-quark prediction. In Fig. ^ on the next page, we plot 
the cross section as a function of Y using the same approximations as in Fig. [6|. Also the 



"in evaluating the analytic high-energy limit, we used the equivalent photon approximation for the lepton 



current (C. 21), as in Ref. [E5f 
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Figure 7: Same as Fig. g, but with the total cross section as a function of Y. 



conclusions are basically the same as for Fig. ||, namely the three curves converge only for 
Y > 7, while for Y < 6 the analytic high-energy limit significantly overestimates the exact 
contribution of the diagrams of the g class. Note that the solid and dashed lines of Fig. ^ 
are the same as the solid and dotted lines of Fig. || but for using the renormalisation scale 
( |4.4| ) instead of ( |4.1| ). Fig. || on the next page has the same content as Fig. 0, but it is 
for a e + e~ future linear collider running at y/s = 500 GeV. For the sake of illustration, we 
have taken the following acceptance cuts: the lepton energies are larger than 40 GeV, the 
lepton tagging angles are between 20 and 70mrad and the hadronic energy is larger than 
20 GeV. On the photon virtualities, the cuts above imply that Q\ 2 > 4 GeV 2 , while the 
average virtualities are (Qf) ~ 36 GeV 2 . About the three curves, the same conclusions as 
for Fig. ^ can be drawn. However, the much larger statistics (at the designed luminosity, 
C = 3.4- 10 34 cm~ 2 s _1 [34], we expect about 1700 events in ten days of continuous running) 
should make also the Y > 7 region available to the analysis. 

In conclusion, we have considered three successive approximations to the total cross 
section at 0(a§): 

• the contribution of the diagrams of the g class only, Eq. (|3.6j); 

• the high-energy limit ( |3.4D of the squared matrix element, integrated over: 
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Figure 8: Same as Fig. 0, but for a future linear collider running at y'i = 500 GeV. 



the exact phase space fl3.8|); 

the high-energy phase space ( 3,14| ); 



and we have seen that, although the contributions of the diagrams with gluon exchange 
in the crossed channel are numerically important in the high-Y or high-IY regions, in the 



kinematic range of the LEP2 experiments the high-energy limit (3.4) is not sufficiently 
accurate. 



5. Phenomenological results 

In Section |||, we have analysed the distributions in rapidity of the final-state partons, 
and their contribution to the total cross section in the large Y or large W regions. As 
expected, we have found that the diagrams of the / class yield a contribution which in 
shape is very similar to the one of the NLO calculation. Then we may argue that, if 
properly counterweighted by the virtual corrections, which at this moment are unknown, 
they would yield a rather minor numerical contribution, since they are an order in a s 
higher than the NLO one. On the contrary, the diagrams of the g class, which are by 
themselves finite and gauge invariant, have a very different shape in Y, becoming more 
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Figure 9: Total cross section as a function of Y, at NLO (dot-dashed line), and at NLO plus the 
0(a 2 s ) contribution of the diagrams of the g class only (solid line). The shaded band has been 
obtained by varying the renormalisation scale from Hq/2 to 2/xq- The points are the experimental 
data from the CERN L3 Collaboration [^]. In computing the error bars, we added their statistical 
and systematic errors in quadrature. 



and more numerically relevant as Y grows. However, we have seen that in the kinematic 
range of the LEP2 experiments the diagrams of the g class must be evaluated exactly, the 
high-energy limit Q3.4| ) being not sufficiently accurate. Thus in this section we shall analyse 
the total cross section as a function of Y through the NLO calculation and/or the diagrams 
of the g class. 

In Fig. ^, we plot the total cross section as a function of Y, at NLO (dot-dashed 
line) and at NLO plus the 0(a 2 s ) contribution of the diagrams of the g class (solid line). 
The shaded band has been obtained by varying the renormalisation scale ( |4.1| ) hq/2 to 
2/zo, with (io given in ( f4.1| ). The points are the experimental data from the CERN L3 
Collaboration [22]. The corresponding numbers are given in Table In computing the 



error bars, we added the statistical and systematic errors of the data in quadrature. We see 
that adding the diagrams of the g class to the NLO calculation decreases the discrepancy 
between data and theory at large Y, however the data still lie above the theory prediction, 
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even allowing for a scale uncertainty on the latter". In addition, we remind that our 
calculation is performed in the massless limit: no mass effect for final-state charm and 
bottom quarks have been included. In Ref. |^9| it was found the masses to decrease the LO 
cross section by 10-15%. A comparable depletion is expected at NLO. In the case of four 
quark production in the analytic high energy limit (defined in the discussion of Fig. 6), the 
masses were found to decrease the cross section by about 20% [27|. This correction should 
provide a lower bound to the exact mass correction in four quark production. Since the 
four quark contribution dominates over the NLO calculation at large Y, and is negligible 
at small Y, we should expect the inclusion of the mass corrections to decrease the solid 
line of Fig. || by about 10-15% at small Y and by at least 20% at large Y. Therefore the 
inclusion of the mass dependence is expected to improve the agreement between data and 
theory at small Y but to widen the discrepancy at large Y. The same considerations apply 
to Fig. [[(], where the total cross section is plotted as a function of W. This was to be 
expected, since in the large Y limit, Y grows linearly with the logarithm of W. 



AW 


L3 data 


NLO 


NLO ± g class 


(GeV) 


da ee /dW (pb/GeV) 


da ee /dW (pb/GeV) 


da ec /dW (pb/GeV) 


5- 10 


0.0747 ± 0.0096 ± 0.0067 


0.0883 _ 00027 


0.0885 _ 00027 


10- 20 


0.0263 ± 0.0024 ± 0.0024 


0300 +0 - 0001 


0305 + - 0003 
u.u3uo _ 00002 


20- 40 


0.0062 ± 0.0007 ± 0.0006 


U - UU0 ' -0.0003 


0064 +u ' UUUfi 
u - uuu -0.0003 


40-100 


0.0014 ± 0.0002 ± 0.0001 


0004 +u uuul 

U - UULF± -0.0000 


0007 +u ' UUU2 
u.uuui _ oooi 


Ay 


L3 data 


NLO 


NLO + g class 




da cc /dY (pb) 


da ee /dY (pb) 


da cc /dY (pb) 


2.0-2.5 


0.315 ±0.048 ±0.028 


n ncR +0.001 
°- 366 -0.001 


0.368™ 


2.5-3.5 


0.184 ±0.018 ±0.017 


2 03 +0 002 
u - zuo -0.001 


0.208™ 


3.5-5.0 


0.085 ± 0.009 ± 0.008 


U - U(U -0.002 


0.0801™ 


5.0-7.0 


0.037 ± 0.006 ± 0.003 


010 +u uul 

u,uiu -0.001 


0-018™ 



Table 1: Differential cross sections in W and Y for the process e + e~ — > hadrons. For the data the 
first uncertainty is statistical and the second systematic. For the theoretical predictions the error 
is given by the renormalization-scale ambiguity. 



6. Conclusions 



In Ref. [29 1, the question had been addressed of whether the LEP2 data for the total 
cross section of 7*7* — > hadrons could be described by a NLO calculation. It was found 
that the NLO analysis described well the data, except at the high end of the hadronic 
energy spectrum. Through the analysis of the inclusive jet and dijet cross sections, different 
kinematic regions were explored, and it was argued that the region of large Y is particularly 
susceptible to large logarithms of type ln(W 2 / /j^) ■ 

" We have also computed the total cross section as a function of Y, by using only the diagrams of the g 
class, i.e. without the NLO contribution, and evolving a s with the one-loop running. For the high-V region 
(Y > 5) the outcome is compatible in shape with the dot-dashed curve, thus it fails as well to describe the 
data. 
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Figure 10: Same as in Fig. || but with the total cross section as a function of the hadronic energy 
W. 



In this work, we have included in the analysis the four parton final states, which are 
part of the O{o 2 s ) contribution to the total cross section. The four-parton final states, 
which have been included in the massless limit, contain the diagrams with gluon exchange 
in the crossed channel, i.e. the diagrams of the g class, which constitute the leading order of 
the BFKL resummation. In Section ||, we have shown that indeed they play an important 
role in the large Y region, however they must be evaluated exactly. In fact, the high-energy 
limit ( |3.4| ), which constitutes the kinematic framework of the BFKL resummation, is not 
sufficiently accurate at LEP2 energies, when compared to the experimental accuracy. Thus, 
if a BFKL resummation is used in the large Y region, we expect the subleading logarithmic 
corrections to be sizeable. 

In Section [|, we have shown that the contribution of the diagrams of the g class to the 
total cross section reduces the discrepancy between the theory and the LEP2 data of the 
L3 Collaboration. However, even allowing for the large scale uncertainty, which is intrinsic 
to the diagrams of the g class since they appear for the first time at 0(a|), the LEP2 data 
still lie above the theory. We remind the reader that in the NLO calculation and in the 
exact four-quark contribution quark mass effects have not been included. The inclusion 
of the mass dependence is expected to improve the agreement between data and theory at 



- 20 - 



small Y but to widen the discrepancy at large Y. Thus, in order to describe accurately 
the data for 7*7* — ► hadrons, mass effects should be included, and eventually in the large 
Y region corrections of an order higher than C(a|) should be considered. 
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Note added in proof 

After the completion of this work, we learned that the ALEPH collaboration has also 
finished its analysis of double tagged events at LEP to measure the hadronic cross section 
in virtual photon-photon scattering |37j. In contradiction to the L3 results, this analysis 
indicates that the NLO QCD prediction is sufficient to describe the data in the high-y 
region, but is unable to predict correctly the measured yields in the low-Y region. One has 
to resolve this apparent contradiction before applying the analysis presented in this paper 
to the ALEPH selection cuts. 

A. Chiral-spinor algebra 

In order to evaluate the production rates, we use helicity amplitudes, defined in terms of 
massless Dirac spinors ip±(p) of fixed helicity, 



1±75 
2 



i/)(p) = 



tp±(p) = (P 



± 



(A.l) 



spinor products 



(pk) = (p~\k+) 



[pk} = (p + \k-), 



(A.2) 



currents 



(i\k\j) = (i 

m+i)\j) = (i 



mr) = (ik) [k 3 \ 
w+f)\n 



(A.3) 



and Mandelstam invariants 



s pk = 2p ■ k = {pk) [kp] , t pkq = (p + k + q) 2 , t pkq i = (p + k + q + I) 2 . (A A) 



B. Four-quark production 



The helicity amplitude featuring two quark pairs and two lepton pairs is 

A(l g ,2 g -;3^%,5^;7 Q ,8 Q ) = 4e 4 «?| Tfcjfa A 8 (l, 2; 3,4, 5, 6; 7, 8) , (B.l) 

for two quark pairs of different flavour, and with {1,2} and {7,8} the quark pairs**, and 
{3,4} and {5,6} the lepton pairs. In the colour-stripped sub-amplitude A$, the fermion 
flavours, momenta and helicities are implicit in the labels. A% is divided into the functions 
ag, bg and c$, 



A$(l q , 2 q ; 3 e , 4g, 5|>, 6ff\ 7q, 8q) = Qf q Qf Q a$(l, 2; 3, 4, 5, 6; 7 
?^ 8 (l,2;3,4,5,6;7,8)+Qj Q 

with Q/ 9(Q) the electric charge of the quark q(Q) of flavour f q (Q) and 



+Q\ Ml, 2; 3, 4, 5, 6; 7, 8) + Q 2 fo c 8 (l, 2; 3, 4, 5, 6; 7, 8) , (B.2) 



o 8 (l, 2; 3, 4, 5, 6; 7, 8) = 5a (l, 2; 3, 4, 5, 6; 7, 8) + <? a (l, 2; 6, 5, 4, 3; 7, 8) 
+ g b (l, 2; 3, 4, 5, 6; 7, 8) + g c (l, 2; 3, 4, 5, 6; 7, 8) 
+ ({l,2}<-{7,8}) 

68(1, 2; 3, 4, 5, 6; 7, 8) = /(l, 2; 3, 4, 5, 6; 7, 8) + /(l, 2; 6, 5, 4, 3; 7, 8) (B.3) 
cb(1, 2; 3, 4, 5, 6; 7, 8) = 6 8 (7, 8; 3, 4, 5, 6; 1, 2) , 

where in the quark pair exchange ({1,2} «-> {7,8}) we swap the momentum and helicity 
labels of the (anti)quarks. In Eq. ( p.2| ) we have factored the flavour dependence in the quark 
electric charges, thus the functions a§, &s an d c§ are independent of the quark flavours. 
In addition, because of the explicit sum in Eq. (|B.3| ) over the different orientations of the 
quark lines, the labels of the partons 1 and 7 refer only to quarks, and not to antiquarks. 
For distinct flavours each of the functions as, b$ and eg is gauge invariant. The functions 
g(f) refer to diagrams which feature gluon (quark) exchange in the crossed channel, Fig. |1[ 
The functions g b and g c are symmetric under the exchange of the quark pairs and of the 
lepton pairs, 

#(1, 2; 3, 4, 5, 6; 7, 8) = 5i (7, 8; 6, 5,4, 3; 1, 2) i = b,c (B.4) 
For the configuration (1~, 2 + ; 3~, 4 + , 5 + , 6~; 7~, 8 + ), the functions g a , g b , g c and / are, 
ga = ■ (13) [58] (6|(5 + 8)|2)(7|(1 + 3)|4) 

•S34 ^56 ^134^568*1234 

gb = - (13) (6 7) [2 8] ([5 6] (6|(1 + 3)|4) + [5 7] (7|(1 + 3)|4)) 

S34 S56 £134 *567 ^1234 

(1 7) [2 4] [5 8] ((2 3) (6|(5 + 8)|2) + (43) (6|(5 + 8)|4)) 
9c = % \o. I ) 

•S34 S56 *234 *568 *1234 



/ 



(13) [25] (6|(2 + 5)|8)(7|(1 + 3)|4) 

S34 S56 ^78 tl34 *256 



*We normalise the colour matrices in the fundamental representation as Xv(T a T b ) — 5 a 
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| (17)[25](6|(2 + 5)|4)(3|(1 + 7)|8) 

■S34 «56 s 78 *178 *256 

| (13)[28](7|(2 + 8)|5)(6|(1 + 3)|4) \ 

■S34 S56 S78 tl34 *278 / 

For all of the other helicity configurations, the functions g a , gi,, g c and / assume a functional 
form which is in principle different, however, the other lepton- helicity configurations are 
simply obtained by exchanging the labels 3 and 4 and/or 5 and 6 in Eq. (|B,l[) . Analogously, 
we show in Appendix BA that the other quark-helicity configurations are obtained by 
exchanging the labels 1 and 2 and/or 7 and 8. 

In the squared amplitude, the sum over distinct flavours can be written as 



2 

JqjQT'fq 



\ l^ 8 (l,2;3,4,5,6;7,8)| 2 

E {Q\Q) Q + 2R <blc 8 )} + 2QlQ fQ Re(a* 8 b 8 ) + 2Q / ,Qj g Re(a$c 8 )} 

fQ¥=fq 

+ ("/-i)(E<#) (n 2 + m 2 ) 



fqjQ¥=fq 

' (B.9) 



/ 

with rif the number of quark flavours, and 

YQ A f = Q^u + Q A d n d (B.10) 

/ 

E QfM Q = QuQ>u(n u - 1) + Q l d Q 3 d n d {n d - 1) + (q^qJ + QiQ^) n u n d 

fqjQ^fq 

with i,j any integer power, and with Q u = 2/3, Q d = —1/3 and n u ^ the number of 
up(down)-type quarks. In Eq. ( |B.9| ) the factor of 1/2 appears in order to avoid double 
counting, since the amplitudes are symmetric with respect to the interchange of the two 
quark lines. 

For two quark pairs of equal flavour, the sub-amplitude A% of Eq. (|B.2j ) becomes 



4 d (l„ 2 ? ; 3 £ , A b h llt 6 r , 7 q , 8 ? ) = Qj4 d (l, 2; 3, 4, 5, 6; 7, 8) , (B.ll) 

with 

4 d (l, 2; 3, 4, 5, 6; 7, 8) = 5a (l, 2; 3, 4, 5, 6; 7, 8) + <? a (l, 2; 6, 5, 4, 3; 7, 8) 

+ g b (l, 2; 3, 4, 5, 6; 7, 8) + g c (l, 2; 3, 4, 5, 6; 7, 8) (B.12) 
+ /(l, 2; 3, 4, 5, 6; 7, 8) + /(l, 2; 6, 5, 4, 3; 7, 8) 
+ ({l,2}<-{7,8}). 

Note that in this instance the diagrams with gluon exchange in the crossed channel, corre- 
sponding to the functions g a , gb and g c , are in the same gauge class, while those featuring 
quark exchange in the crossed channel form a different gauge class. In addition, quarks of 



- 23 - 



equal flavour are indistinguishable, thus we must add to Eq. (|B.1| ) the contribution with 
the quarks (but not the anti-quarks) exchanged, and antisymmetrise the whole amplitude 
in the colour and momentum labels, 



At(l q , 2 q ; 3 £ , 4 b 5e, 6i>; 7 q ,8g) = Ae*g 2 s 



x 2 [( P )n« 2 rU + ( r )« ? « 2 rkiJ [4 d (l,2;3,4,5,6;7,8) - <(7, 2; 3, 4, 5, 6; 1, 8) 



+ [(T a ) lll2 (T a ) tTh - (T a ) lrh (T a ) tlh ] [Af(l, 2; 3, 4, 5, 6; 7, 8) + A$(7, 2; 3, 4, 5, 6; 1, 
4e 4 5l [(r a ) n 5 2 (T a )^ 8 4 d (l,2;3,4,5,6;7,8) - (T%; 2 (^^(7, 2; 3, 4, 5, 6; 1, 1 



(B.13) 



i.e. we can antisymmetrise Eqs. ( [B.lD and ( |B.ll ) by subtracting the same expression with 
the colour and momentum labels of the quarks exchanged. 

In crossing to the physical region, we choose 4 as the incoming electron and 6 as 
the incoming positron. For a fixed lepton-helicity configuration, e.g. (3^,4^,5^,6^), the 
production rate is obtained by summing over the quark-helicity configurations, 

da(3J,4f, 5 J, 6^) = ^dV 6 16 (N 2 C - l) (47ra em ) 4 (47ra s ) 2 



,2+;3-,4 + ,5 + ,6-;7- 



+ |^ 8 (2-,l+;3-,4+,5 + ,6-;7- 



+ |A 8 (l-,2+;3-,4 + ,5 + ,6-;8-,7 + )| 2 + |yl 8 (2-,l + ;3-,4 + ,5 + ,6-;8-,7+)| 2 
+ \{T,Q 4 f) [l4 d (l" ) 2+;3-,4+,5 + ,6-;7-,8 + )| 2 + |4 d (7^,2+;3-,4+,5+,6-;l 



+— Re[4 d (l-,2+;3-,4 + ,5 + ,6-;7-,8 + )X d (7-,2+;3-,4 + ,5 + ,6-;l- 



+l4 d (2-,i H 

+|4 d (7",2^ 

+l4 d (2-,i H 



3-,4 + ,5+,6-;7- 



+ |4 d (l-,2 + ;3-,4 + ,5 + ,6 



3-,4 + ,5+,6-;8-,l + )| 2 + |4 d (2-,7 + ;3-,4 + ,5 + ,6- 
3-,4+,5 + ,6-;8-,7 + )| 2 + |4 d (2-,7 + ;3-,4 + ,5 + ,6- 



1- 



+-^Re [4 d (2~, 1+; 3-, 4 + , 5 + , 6-; 8-, 7+)*a 8 d (2-, 7 + ; 3-, 4 + , 5+, 6-; 8-, 1+)] 



* „id/ 



N, 



where dV% is the phase space for the e + e qqQQ final state, 



(B.14) 



(B.15) 



with i = 1, 2, 3, 5, 7, 8. In Eq. ( B.14 ) we have performed explicitly the sum over quarks of 
equal flavour, and we have multiplied by the symmetry factor 1/4 for two identical quarks 
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and two identical antiquaries. The sum over quarks of different flavour is performed in 
Eq. ( p3.9|) . The unpolarised rate is given by averaging over the rates for the four lepton- 
helicity configurations. 

B.l Symmetries under helicity flips of the quark lines 

Symmetry relations between the functions g a , gb, gc and / with respect to the configuration 
(1- 2+;3- 4+,5+ 6-;7-,8+), 

(a) under helicity flip of the pair {1,2}, 



(B.16) 
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(b) under helicity flip of the pair {7, 8}, 



9a(l~ 


2+ 


3" 


,4 + ,5+,6" 


7+ 


,8" 


") = -06(1- 


,2+ 


;3" 


,4+,5+,6- 


;8-,7+) 


9a(7 + 


8" 


3" 


,4+,5 + ,6" 


1" 


,2 4 


") = -9c(l- 
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,5+,4+,3- 


;8-,7 + ) 


9b(l- 


2+ 


3" 


,4 + ,5 + ,6" 


7+ 


,8' 


') = -9a(l~ 
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,4+,5 + ,6" 


;8~, 7+) 


9c(l- 
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3" 


,4+,5+,6~ 
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,8" 
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8", 7+) 
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8~ 


3~ 


,4+,5+,6~ 


1" 


,2 4 


= ~/(8- 


7+ 
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5+,4 + ,3"; 


1-.2+) 



(c) under helicity flips of the pairs {1,2} and {7,8}, reflection 



9a(l + 


2" 


3" 


,4 + 


5+,6~ 


7+,8~ 


) = 5a (8- 


-,7+;6- 


-,5+4+3- 


;2~,i + ) 


9a(7 + 


8" 
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(B.17) 



(B.18) 



Eqs. (|B.16|) - (|B.18| ) show that the helicity flip of the pairs {1,2} and/or {7, 8} reshuffles 
the functions in Eq. (^^) however it does not change their sum but only the overall sign, 
thus the helicity flip of the pairs {1,2} and/or {7,8} of the sub-amplitude (|3.3|) can be 
achieved by merely exchanging the corresponding labels. 



C. The impact factor for j*g* —> qq 

In order to derive the impact factor for 7*5* 
amplitudes for e + e~ — * qqgg or those for e+e" - 



-» qq, the simplest is to use either the 
qqQQ. They are collected, for instance, 
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in Ref. [35]. We shall take a lepton and a parton (quark or gluon), of momenta p a and pb 
respectively as the incoming particles, and a lepton, a quark pair and a parton of momenta 
Pa 1 : Pi: P2 and as the outgoing particles. All momenta are taken as outgoing, so that 
momentum conservation reads Pa+Pb+Pa' +P1+P2+P0 1 = 0. We use light-cone coordinates: 
p± _ p0 ±pz anc [ com pl e x transverse coordinates p± = p x + ip y , and write a four vector 
= (p ,^) 1 ,^ 2 ,^ 3 ) as (p + ,p~',p±,p±), where p± = p x — ip y . With this notation the scalar 
product is 2p ■ q = p + q~ + p~q + — p±q± — p±q±- For real four vectors, e.g. for momenta, 
for which p±_ = Pj_, we shall use the shorter notation p = {p + ,p~ ;p±). 

We define the photon momentum to be k = —p a — Pa' , and choose as the reference 
frame the virtual photon-parton frame, defined by pb = (0,p^;0^) and k = (k + , k~; 0j_). 
Momentum conservation requires that 

k + = p + + p + + p+ 

k~ = Pi +P2 +Pb + Pb' (C- 1 ) 
= Pi± + P2± + Pb'L 

Now we take the high-energy limit, where the outgoing partons are strongly ordered on 
the light cone and have comparable transverse momentum, 

pt-Pt^Pb>> Pi-P2<-Pv'i IjPulI - |P2±| - |P6'±| • (C2) 



While the transverse components of Eq. ( |C . 1| ) remain untouched, the light-cone components 
are approximated by k + ~ p^ + p\ and k~ ~ p^ +p b ~,. In addition, the virtual photon- 
parton centre-of-mass energy 



(k - p b ) 2 = k + k~ - k+pl ~ k + p y (C.3) 



is required in the high-energy limit to be much larger than the virtual photon momentum 
transfer k 2 , i.e. s 7 * p 3> \k 2 \. This entails that p^ 3> \k~\. Thus in the momentum 
conservation along the minus direction the momentum k~ can be neglected, and we can 
summarise the momentum conservation in the high-energy limit as 

-pi - pi = k+ - pi +pt 

-Pb - Pb' (°- 4 ) 

= Pl± + P2± + Pb'± ■ 

Eq. ( |C.4D can be viewed as defining two scattering centres through the + and — momen- 
tum conservation, which act independently. The two scattering centres are linked by the 
transverse momentum conservation only. 

Next, we must approximate the exact amplitudes in the high-energy limit. The colour 
decomposition of the amplitude for e + e~ — ► qqgg in the conventions of Ref. [35], is 

Mh,b,b',2 q ;a E ,a' e ) = -2e 2 g 2 s Q fq £ {T a °(»T a ^')) ilh A 6 (l q , cr{b), a(b'), 2 ? ; a 8 , a! e ) 

creS 2 

(C.5) 
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with sub-amplitudes^ 



A 6 (l + ,b + ,b' ,2-, a', a 



1 



(1/2) [2 6] [la'] (a\{b' + 2)\b) 
[V 2] t W2 



(b'l) [16] (2a) (b'\(l + b)\a') 

(a\(b' + 2)\b)(b>\(l + b)\a>) 
(lb) [b'2] 



(C.6) 



,4 6 (l + ,6~,&' + ,2",a~,a' 



Sbb'Sa 



lb'] 2 (2 a) <6|(l + 6')|a') 



[16]* 



166' 



(C.7) 



+ (62)' [la'] (a|(6 + 2)|6') + [1 6'] (6 2) [1 a'] (2 a) 



y2>tw 2 [16] (V 2) 

The sub-amplitudes ( |C.6| ) and ( p,7|) are symmetric under the exchange 

1 <-» 2 , 6 6 , a <-> a' , (r/) <-> [ji] . 

Alternatively, the amplitude for e + e _ — ► qqQQ can be used, 

8 

^ 6 ree (l g ,6 ,6^2 ? -;ae,a' e ) = -2e 2 5 | Y, T hn T iy% 

a=l 

x [<2/ 9 ^6 ree (l? ; 6q, 6'q, 2 5 ; a g , a' e ) + Q/ Q ^^(b'Q, 2 ? , 1 9 , 6q; o g) a' e 

with sub-amplitude 

4 ree (l+,6+,6'-,2-a-,a' + ) 

16] (a 2) (6'|(1 + 6)|a') (6' 2) [a' 1] (a|(6' + 2)|6) 



(C.8) 



(C.9) 



Sbb'Sa 



hbb' 



+ 



hb'2 



(CIO) 



In order to evaluate the sub-amplitudes above, we need to compute the Mandelstam in- 
variants and the spinor products in the high-energy limit. In the virtual photon-parton 
frame, the three-particle invariants can be written exactly as 



Uaa> = (Pi + Pa + Pa' 



k 2 



^ + I 1 2 Pi h.2 
Pi k 



* = 1,2, 



(C.ll) 



where = pf /k + , with i = 1, 2, are the momentum fractions of the final-state quarks with 
respect to the virtual photon. In the high-energy limit, we can fix x\ = x and x 2 = 1 — x, 
and rewrite the invariants ( |C.11| ) as 



haa' = hbb' = (1 - x)k 



2 \Pl±\ 



haa' = hbb' = Xk 



\P2±\ 

(1-x) 



(C.12) 



^We use the same sub-amplitudes as in Ref. |35|, but we ignore the overall factor of i. In addition, we 
have neglected the configuration with like-helicity gluons, since it subleading in the high-energy limit. 
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In light-cone coordinates, a generic spinor product can be written as 



P 



I Pi 



(PiPj) = Pi±\Mr ~Pj±\ Mr 



Pi 



(C.13) 



In the kinematics ( C.4 ), the spinor products are 

(PaPb) = 

(PaPb'} = 



(Pt + V2 + Va')Pv ' 



py 

~ i^-^(PaPb) , 




(C.14) 



(PaPj) 

(PjPb) 
(Pb'Pb) 

(PjPb'} 



with j = 1,2, a' , and where we have taken the phase conventions of Ref. [36]. Momentum 
conservation ( C.4 ) implies that p£ = —p^, — pf — p\ and p a ± = —p a '±- 

Using the invariants ( |C.12[ ) and the spinor products ( |C,14j ) , the high-energy expansion 
of Eq. ( |C.5| ) can be written as 



A 6 (l»,br,b' p ',2g\a^\a> X e ) 



2s 



TP 



9se 2 Q fq If A V2V f (a^; a' e \ 1J, 2? 



t 



-^g s f bb ' c C(b^) 



, (C15) 



where s 7 * p ~ k + p^ is the virtual photon-parton centre-of-mass energy, t = sw — — \Pb'±\ 2 = 
— \q_\_\ 2 is the momentum transfer, with q = Pb+Pb' the momentum of the gluon exchanged 
in the crossed channel, and where A , v and p denote the helicity of the lepton pair, the 
quark pair and the gluons, respectively. Using the impact factor for the gluon vertex, 
C(b^;b'~) = Pb'±/pl'±, according to the conventions of Ref. fl36|P we obtain the impact 
factor for eg* — > eqq, 



Vf(a g ;a' e , l+,2 



(C.16) 



i 



+ 



x(l — x) 



-x(l — x)k 2 + |j»2±| 2 

x(l — x) 
-x(l — x)k 2 + |pij_| 2 



Pa± 



Pl± 




+ Pa> 



Pl±\ 




-pi 



Pi 



H In Ref. |^] the generators of the group are normalised to 1/2, while in this paper to 1. We introduced 
the explicit l/v2 factor in Eq. (C.15) to take into account this difference. 
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Note that in Eq. ( |C.16| ) there is no divergence when the momenta p\ and P2 of the quark 
pair become collinear. As a check of the calculation (or of high-energy factorisation), we 
have evaluated also the amplitude ( p.9|) in the high-energy limit, and have obtained the 
same result as in Eq. ( C.15| ), up to the substitution 



ig s f bb ' c C{b»» ; &7 ) <-> g s T b % C(bp' ; b> u f ) , (C. 17) 

with impact factor for the quark vertex, C(b^;b'~) = (— i)pb'±/\Pb'±\- 

Using discrete symmetries of the helicity amplitudes, we find the impact factors with 
the helicities of the quark and/or the lepton pairs flipped, 

VKa±;a'- 1-2+) = [V f (aJ;a't, 1+ 27)]* , 

V f (a=-;a / + 17,2+) = V f (aJ; a' + , 2+ 17) , (C.18) 



F / (a+;a / 7;l+ 27) = [^(aj; a'+; 2+ 17)]* . 



The impact factor for the backward kinematics is 

V b (a^;a't; 1+ 27) = e i( ^ ^/(a'g ;a+; 2+ 17)|+^_ , (C.19) 

where the phase e 1 ^ for us is immaterial since we compute squared amplitudes and the 
index + — > — means that the plus components in Eq. ( p,16| ) are replaced with minus 
components. 

Eq. ( |C16 ) can be decomposed further in terms of an impact factor for j*g* — > qq by 



factoring the lepton current times the photon propagator {p a \il^\Pa'){—'i'9nv)/k\ 

F / (a7;a'+ 1+ 2r) = ( a - \yf\ a '-)^-V$ (k; 1+ 27) . (C.20) 
In light-cone notation, the lepton current is 



(Pa ~ \i-f\p*-) = "2 J-ptpt„ -mA, Pl'J^ , (C21) 



-T&T& X Pa V Pa' 

Pa P a i ' ' a 

with the 7 matrices chosen in the chiral representation as in Ref. [36|. The impact factor 
Vj ) i(k; 1 + , 27) for the 7*3* — ► process can be written as 

V£(k;l+, 27) (C.22) 



-x(l — x)/c 2 + |^2±| 2 



(x(l - x)k + , k ; xp 2 ±, -(1 - x)p* 2 ^j 



y/x(l-x) ( + \pi±\ 2 h - n x * 

' x(l - x)k^, — rj-k ; -xpu_, (1 - x)p 1± 



—x(l — x)k 2 + \pi±\ 2 \ ' — k 2 

Using k = (k + , 0j_), we can easily check that k^Vji{k\ 1 + , 27) = 
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Finally, if we contract with the polarization vector, = {e + , e ; e±), of the virtual 
photon, we obtain 



V^(k;l+,2- q ) 



(C.23) 



-x(l — x)k 2 + \p2±\ 



a/x(1 — x)k + e + 



\P2±l 

—k 2 



-k e + + (1 - x)p* 21 _£± - xp 2 j_e* L 



^xJT 



\j x(\ — x)k^ 



£ + 



\PlM 



-k e + - (1 - x)p* 1± e± + xpi^s^ 



—x(l — x)k 2 + \pi±\ 2 \ —k 2 
We have checked that using this form of the impact factor, we can reproduce the cross 



section in the high-energy limit obtained in Ref. [25]. 

The impact factor with the helicity of the quark pair flipped is obtained by exchanging 
the momenta of the quark and antiquark, 



<.(*;!- 2+) = 2+ K 



(C.24) 



The impact factors for the backward kinematics, k~ 3> k + , have the same functional form 
as given by Eq. flC.22j ) (up to a phase) with the + and — components as well as the quark 
helicities interchanged. 

As a final check, we computed the impact factor for eg* — ► eqq, Eq. ( C.16| ), also in the 
electron-parton frame, and verified that it agrees with the colour-subleading piece, termed 
-E?2, of the impact factor for qg* — ► qQQ, computed in Ref. [36]. 
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